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Abstract 

The exclusive electroproduction of J/ip mesons, ep — > ep J /if), has been studied 
with the ZEUS detector at HERA for virtualities of the exchanged photon in 
the ranges 0.15 < Q 2 < 0.8 GeV 2 and 2 < Q 2 < 100 GeV 2 using integrated lumi- 
nosities of 69 ptT 1 and 83 pb" 1 , respectively. The photon-proton centre-of-mass 
energy was in the range 30 < W < 220 GeV and the squared four-momentum 
transfer at the proton vertex \t\ < 1 GeV 2 . The cross sections and decay an- 
gular distributions are presented as functions of Q 2 , W and t. The effective 
parameters of the Pomeron trajectory are in agreement with those found in J/ip 
photoproduction. The spin-density matrix elements, calculated from the decay 
angular distributions, are consistent with the hypothesis of s-channel helicity 
conservation. The ratio of the longitudinal to transverse cross sections, ctl/&t, 
grows with Q 2 , whilst no dependence on W or t is observed. The results are in 
agreement with perturbative QCD calculations and exhibit a strong sensitivity 
to the gluon distribution in the proton. 
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1 Introduction 



The exclusive electroproduction of light (p,cu,<f>) and heavy (J/tf},i/;', T) vector mesons, 
ep — > eVp, has been investigated at HERA [1-9]. The increased precision of the recent 
data allows the study of the dependence of this process on the different scales involved: 
the mass squared of the vector meson, My, the square of the centre-of-mass energy of the 
photon-proton system, W 2 , the exchanged-photon virtuality, Q 2 , and the four-momentum 
transfer squared at the proton vertex, t. 

Exclusive electroproduction of vector mesons involving a sufficiently large scale is calcula- 
ble perturbatively because of the QCD factorisation theorem [10]. QCD-based models of 
this process assume that the exchanged virtual photon, seen from the proton rest-frame, 
fluctuates into a quark-antiquark pair which interacts with the proton via the exchange 
of two gluons in a colour-singlet configuration. After the interaction, the qq pair becomes 
a bound state. The cross section is proportional to the generalised parton distribution 
functions (GPDs) [11] of the proton, which contain information on the momentum distri- 
butions of the partons in the proton and their correlations. At the leading-order approxi- 
mation in ln(l/x) and vanishing t, the generalised gluon distribution can be approximated 
by the usual gluon distribution. The gluon density is probed at x ~ (Q 2 + M v )/W 2 and 
at a scale p 2 ~ Q 2 + M v [12]. The cross section is thus expected to rise steeply with W, 
a reflection of the steep rise of the gluon density as x decreases. 

Data from exclusive p production [1-3] show that the cross section u{^f*p — > pp) rises 
with W as W 5 , with 5 increasing with Q 2 from about 0.2 at Q 2 = (photoproduction) to 
about 0.8 at Q 2 « 30GeV 2 . However, in the case of exclusive J/ip production the cross 
section rises steeply with W even for photoproduction [9]. It is therefore interesting to 
investigate J/i[) production at larger values of Q 2 . 

This paper presents measurements of the exclusive electroproduction of J/ip mesons. 
Cross sections are given as functions of W, Q 2 and t. The W dependence is also studied 
as a function of t. The helicity structure of the J/i/j has been investigated to test s- 
channel helicity conservation (SCHC) and to extract the ratio of the cross sections for 
longitudinally (cr^,) and transversely (err) polarised virtual photons, R = cfl/^Ti as a 
function of W, Q 2 and t. The results are compared to perturbative QCD (pQCD) model 
calculations. 

The data cover the kinematic range 30 < W < 220 GeV and |t| < 1 GeV 2 for two ranges of 
photon virtuality: 0.15 < Q 2 < 0.8 GeV 2 (low-Q 2 sample) and 2 < Q 2 < 100 GeV 2 (high- 
Q 2 sample). The low-Q 2 sample was measured in the e + e~ decay channel and the high-Q 2 
sample in both e + e~ and p + p~ channels. The low-Q 2 range has been measured for the 
first time. The high-Q 2 sample represents more than an order of magnitude increase in 
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statistics compared to the previous ZEUS results [1], and extends both the W and Q 2 
ranges of the measurement. 

2 Experimental set-up 

The data used for this measurement were taken at the HERA ep collider using the ZEUS 
detector in 1998-2000. During this period, HERA operated with a proton energy of 
920 GeV and an electron 1 energy of 27.5 GeV. The data correspond to integrated lumi- 
nosities of 69 pb -1 for the low-Q 2 sample and 83 pb" 1 for the high-Q 2 sample. 

A detailed description of the ZEUS detector can be found elsewhere [13]. A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles were reconstructed in the central tracking detector (CTD) [14] covering 
the polar-angle 2 region 15° < 9 < 164°. The transverse-momentum resolution for full- 
length tracks is a(p T )/pT = 0.0058p T © 0.0065 © 0.0014/p T , with p T in GeV. 

The high-resolution uranium-calorimeter (CAL) [15] consists of three parts: the forward 
(FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part is subdivided 
transversely into towers and longitudinally into an electromagnetic section (EMC) and 
either one (RCAL) or two (FCAL and BCAL) hadronic sections. The CAL covers 99.7% 
of the total solid angle. The energy resolution obtained from test-beam measurements 
was a(E)/E = 0.18/y/E in the electromagnetic sections and a(E)/E = 0.35/^/E in the 
hadronic sections, with E in GeV. 

The forward plug calorimeter (FPC) [16] was a lead-scintillator sandwich calorimeter with 
readout via wavelength shifter fibres. It was installed in the beamhole of the FCAL and 
extended the pseudorapidity coverage of the forward calorimeter from rj < 4 to rj < 5. 

The beampipe calorimeter (BPC) [17] was a tungsten-scintillator sampling calorimeter 
installed to measure scattered electrons at small angles, in the range 1.15° < 180° — 9 < 
2.30°. The energy resolution, as measured under test-beam conditions, was a(E)/E = 
0.17 /\/E, with E in GeV. The impact position of the scattered electron was measured 
with an accuracy of about 0.5 mm. 

The small-angle rear tracking detector (SRTD) [18] consists of two planes of scintillator 
strips read out via optical fibres and photomultiplier tubes. It is attached to the front 

1 Hereafter, both e + and are referred to as electrons, unless explicitly stated otherwise. 

2 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
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face of the RCAL and covers an angular range between 4° and 18° around the beampipe. 
The SRTD provides a transverse position resolution of 0.3 cm [9]. 

The hadron-electron separator installed in the RCAL (RHES) [13] consists of silicon diodes 
placed at a longitudinal depth of three radiation lengths. The RHES provides an electron 
position resolution of 0.9 cm for a single hit and 0.5 cm if the shower spans at least two 
adjacent pads [19]. 

The luminosity was determined from the rate of the bremsstrahlung process ep — > e'-fp, 
where the photon was measured with a lead-scintillator calorimeter [20] located at Z = 
-107 m in the HERA tunnel. 



3 Kinematics and cross sections 

The following kinematic variables are used to describe exclusive J / i/j production, e(k)p(P) - 
e(k') J/i(j(v)p(P'), where k, k', P, P' and v are, respectively, the four-momenta of the in- 
cident electron, scattered electron, incident proton, scattered proton and J/ip: 

• Q 2 = —q 2 = —(k — k') 2 , the negative four-momentum squared of the virtual photon; 

• W 2 = (q + P) 2 , the squared invariant mass of the photon-proton system; 

• y = (P ■ q)/(P ■ k), the fraction of the electron energy transferred to the proton in the 
proton rest frame; 

• x = Q 2 /(2P ■ q), the Bjorken variable; 

• t = (P — P') 2 ., the squared four-momentum transfer at the proton vertex. 

The kinematic variables were reconstructed with the "constrained" method [1] which uses 
the momentum of the J/ ip and the polar and azimuthal angles of the scattered electron. 

The ep cross section can be expressed in terms of the transverse, cr^, and longitudinal, 
o"l, virtual photoproduction cross sections as 

d 2 a ep ^ e J l^ v 

dydQ 2 =^(y,Q 2 )(a T + ea L ) 1 

where IV is the flux of transverse virtual photons [21] and e is the ratio of longitudinal 
and transverse virtual photon fluxes, given by e = 2(1 — y)/(l + (1 — y) 2 )- In the kinematic 
range studied here, e is in the range 0.8 < e < 1, with an average value of 0.99. 

The virtual photon-proton cross section, a 1 * p ^ J ^' p = or + ecr^, can be used to evaluate 
the total exclusive cross section, <J^ Q ^ J ^ P = <j t + a L , through the relation 

■y*p->J/i>p _ _ I" R „j*p-+J/ii>p 
~ 1 + eR 
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where R = olI°t is the ratio of the cross sections for longitudinal and transverse photons. 
The helicity structure of the J/ip production is used to determine R as described in 
Section mm 

4 Reconstruction and selection of the events 

The signature of exclusive J /if) electroproduction, ep — > e J/tjjp , consists of the scattered 
electron and two charged leptons from the J/ip decay, e + e~ or The scattered 

proton is deflected through a small angle and escapes undetected down the beampipe. 

The events were selected online by a three-level trigger [22,23]. For the low-Q 2 sample, 
the trigger [9] for J/ip photoproduction events with decay to the e + e~ final state was 
used, while for the high-Q 2 sample, the trigger required a scattered electron in the CAL 
with energy greater than 4 GeV. 

The following criteria were applied offline to reconstruct and select the events [24]: 

• the energy and position of the scattered electron were measured in the BPC for the 
low-Q 2 sample and in the CAL for the high-Q 2 sample. The energy was required 
to satisfy E e > 10 GeV. The position measurement of the CAL was improved using 
the SRTD (88% of the events) and the RHES (10% of the events). To ensure full 
containment of the electromagnetic shower, fiducial cuts were applied to the impact 
position of the electron on face of the RCAL; 

• the J/ip mesons were reconstructed from the decay leptons. Two tracks of opposite 
charge, well-reconstructed in the CTD with > 0.2 GeV, were selected (two-track 
events). In the case of the electron decay channel, events were also selected by requiring 
one well-reconstructed CTD track and one CAL energy cluster [25] not related to the 
track (one-track events). In addition: 

— the two-track events were required to have the higher-momentum track matched 
to a calorimeter energy cluster for which the fraction of the energy deposited in 
the EMC was consistent with that of an electron or a muon; 

— the one-track events were accepted if, in addition to the measured CTD track as- 
sociated with a CAL cluster, the second cluster lay in the angular range outside 
the CTD acceptance with energy between 2 and 10 GeV. Both clusters were re- 
quired to have a fraction of energy deposited in the EMC consistent with that of 
an electron. 

For both types of events, one additional CTD track was allowed. If present, this track 
was required to match the scattered electron. Events with further tracks were rejected; 
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• the position of the reconstructed vertex was required to be compatible with that of an 
ep collision; 

• to remove events with large initial-state radiation and to reduce the background from 
photoproduction, the requirement 45 < 5 < 65 GeV was imposed, where 5 = ^^(1 — 
cos^j), Ei is the energy of the i th calorimeter cell, 0j its polar angle and the sum runs 
over the cells associated to the final-state leptons; 

• to suppress non-exclusive events, the energy of each CAL cluster not associated to any 
of the final-state leptons was required to be less than 0.3 or 0.4 GeV, depending on 
the CAL section; these thresholds were set to be above the noise level of the CAL. To 
suppress further the contamination from proton-dissociative events, ep — * eJ/tpY, the 
energy in the FPC was required to be less than 1 GeV and the sum of the energy in 
the FCAL cells surrounding the beamhole to be less than 0.5 GeV. These cuts restrict 
the mass of the proton-dissociated system, Y, to My < 3.0 GeV. 

Unless otherwise stated, the results are quoted in the following kinematic range: \t\ < 

1 GeV 2 , 30 < W < 220 GeV for the electron channel and 45 < W < 160 GeV for the muon 
channel. The larger W range for the electron channel was achieved by the inclusion of 
the one-track events. The Q 2 range was 0.15 < Q 2 < 0.8 GeV 2 for the low-Q 2 sample and 

2 < Q 2 < 100 GeV 2 for the high-Q 2 sample. 

The final high-Q 2 sample contains 728 events in the muon channel and 955 events in the 
electron channel, 275 of which are reconstructed using only one track. The final low-Q 2 
sample contains 137 events in the electron channel, 16 of which are reconstructed using 
only one track. The distribution of the events in the x-Q 2 plane is shown in Fig. ^ 

5 Monte Carlo simulation 

The acceptance and the effects of the detector response were determined using samples of 
Monte Carlo (MC) events. All generated events were passed through the standard ZEUS 
detector simulation, based on the Geant 3.13 programme [26], and the ZEUS trigger 
simulation package. 

The exclusive process ep — » e J ftp p was modelled using the Zeusvm [27] MC generator 
interfaced to Heracles 4.6.1 [28] to account for first-order QED radiative effects. The 
7* p — >• J/ip p cross section was parametrised as W s e~ b ^(M 2 ^ + Q 2 )~ n . The parameter 
values n = 2.5(2.3), 5 = 0.75(0.7) and b = 4.5(4.5) GeV -2 were used to describe the 
high-Q 2 (low-Q 2 ) data. The leptonic decay of the J ftp was simulated by the Photos 
programme [29] which includes final-state radiation from the decay leptons. This genera- 
tor assumes SCHC and that the ratio of the cross sections for longitudinal and transverse 
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photons is R = 0.5 • (Q 2 /M 2 ^). 

Proton-dissociative events, ep — > eJ/tpY, were modelled using the generator Epsoft [30]. 
The 7* p — > J/ip Y cross section was parametrised as 

dt dM2 oc W*e-^(M^ + Q r n M^ (1) 

with the parameters n = 2.5, S = 0.75, b = 0.81 GeV -2 and (3 = 2.57 chosen as described 
in Section r 



The QED background stemming from two-photon lepton-pair production 7*7 — > l + l~, 
where the virtual photon originates from the electron vertex and the second photon is 
radiated off the proton, was simulated using the Lpair [31] generator at low Q 2 and the 
Grape-Dilepton 1.1 [32] generator at high Q 2 . The QED-Compton-like processes with 
internal photon conversion were also generated with Grape. 



6 Extraction of the J/ip signal 

Figure El shows the invariant-mass distributions of the muon and electron pairs, obtained 
after the selection described in Section 0] The MC distributions for exclusive J/i/j produc- 
tion and the QED background are also shown. The width of the resonance is dominated 
by the detector resolution, which deteriorates at low and high values of W. 



6.1 Non-resonant background 

The non-resonant background was estimated from the MC distributions of the QED- 
background processes: two-photon lepton-pair production and Compton scattering. For 
the low-Q 2 sample, the normalisation of the QED-background was estimated from a two- 
parameter fit of the signal and the background MC distributions to the invariant mass 
spectra of the data. For the high-Q 2 sample, the normalisation was based on the known 
cross sections and the integrated luminosity of the data. After subtraction of the non- 
resonant distributions, the J / if) signal was determined by counting the events in the mass 
windows 2.8 < M M + M - < 3.4 GeV for the muon channel and 2.6 < M e + e - < 3.4 GeV for 
the electron channel. The lower limit on M e + e - was chosen to include events with reduced 
invariant mass due to bremsstrahlung. The contribution of the non-resonant background 
in the signal range is typically 22% for the electron channel and 14% for the muon channel. 

For the high-Q 2 sample, additional background from pions misidentified as electrons or 
muons was studied using a sample of events with two tracks, neither of which were iden- 
tified as a muon or an electron. The contribution was (2.7 ± 0.6)% for the electron 
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channel and (0.8 ± 0.3)% for the muon channel and was subtracted bin-by-bin for the t 
and decay-angle distributions and globally for the W and Q 2 distributions. 

Events from exclusive ip(2S) production contribute to the J/ip sample through two dif- 
ferent decay channels: (i) ip(2S) — > J /tp + neutrals (branching ratio (23.9±1.2)% [33]), 
where the neutrals are not detected in the CAL, and (ii) ip{2S) — > (branching ratios 
(7.3 ± 0.4) • 10" 3 for the electron channel and (7.0 ± 0.9) • 10" 3 for the muon decay chan- 
nel [33]), because of the limited resolution in the reconstruction of the invariant mass. 
The contribution from both these processes to the J / ip sample was determined using MC 
samples under the assumption that a{jp{2S)) / a(J / ip) = 0.166±0.013 [8]. A contribution 
of (1.8±0.2)% was subtracted. 

6.2 Proton-dissociative background 

The remaining source of background consists of J / i/j production accompanied by proton 
dissociation, ep — ► e J/ip Y, where the particles from the breakup of the proton are not 
detected. 

Proton-dissociative events were studied using a sample of diffractive events selected as 
described in Section with the following exceptions: 

• the elasticity criterion (last criterion in Section H} was not applied to the FPC and to 
a region of FCAL of approximately 50 cm radius around the beampipe; 

• events with decay-lepton tracks at angles smaller than 30° with respect to the outgoing 
proton direction were removed to ensure a rapidity gap between the J /if) and the 
system Y. 

Proton-dissociative events were selected by requiring an energy larger than 1 GeV in the 
FPC. The sample of data tagged by the FPC contained 100 events for \t\ < 3 GeV 2 in 
the kinematic range 45 < W < 160 GeV and Q 2 > 2 GeV 2 . The parameters (see Eq.©) 
that best describe the Q 2 , W and t dependences are n = 2.57 ± 0.09, 5 = 0.61 ± 0.40 
and b = 0.81 ± 0.25 GeV -2 . The MC distribution of M Y was tuned to describe the FPC 
energy distribution, yielding f3 = 2.57±0.67. The values for n and 5 are in agreement with 
those described in Section 17.31 The latter are more precise and were used in Epsoft. 
The values of b and (3 are in agreement with those found for proton-dissociative J ftp 
photoproduction [9]. 

The fraction of proton-dissociative events in the elastic sample, / p _di S s, was determined 
from the relation / p _di SS = /fpcCV 6 ' — 1)> where fppQ denotes the fraction of the proton- 
dissociative sample tagged by the FPC and e' = 32% is the FPC tagging efficiency, 
estimated using Epsoft. The fraction of proton-dissociative events in the final sample, 
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averaged over t for |t| < 1 GeV 2 , was / p _diss = (14.2±2.0(stat.)lf;|(syst.))%, independent 
of W and Q 2 . The systematic uncertainty was dominated by the uncertainty on (3. The 
fraction increases from 4% for < |t| < 0.1 GeV 2 to 20% for 0.2 < \t\ < 1.0 GeV 2 . The 
cross sections presented in the next sessions were corrected for this background in bins of 
t, and globally in W and Q 2 . 

7 Results 

7.1 Measurement of cross sections 

In each bin of a kinematic variable, the ep cross section was extracted for each decay 
channel using the formula 

ep-^e J/ij) p (^data ~ -Vbgd)(l ~ Jp-diss) 

ABC 

where iVdata is the number of events in the data and iVbgd is the number of events from the 
non- resonant background (QED processes and pionic background) and i/)(2S) production. 
The overall acceptance is denoted as A, B accounts for the J/ip decay branching ratios [33], 
(5.93 ± 0.10)% and (5.88 ± 0.10)% for the electron and muon channels, respectively, and 
£ is the integrated luminosity. 

The total exclusive photon-proton cross section was calculated as 

atr J/4 ' P = (imQ 2 ,W))d 2 a ep - ieJ ^ p /dQ 2 dW, 

where the effective photon flux $ [34] contains the corrections for bin-centring and R, 
both estimated from the MC simulation. The final cross section was the error-weighted 
average of the cross sections for each decay channel. 

The cross sections are quoted at the QED-Born level. The radiative corrections range 
from 1% to 10% (on average 5%), depending on the kinematic region. 

The cross sections were measured for \t\ < 1 GeV 2 . Assuming da/dt oc e _fe '*', with b = 
4.5 GeV -2 , the correction factor needed to extrapolate to the cross section integrated 
over the full t range is 1.012. In addition, for x > 0.01 both the acceptance and the 
expected cross section are small, and the measurement in this region therefore involves 
an extrapolation, made in order to quote the measurement in bins of W and Q 2 . The 
uncertainty introduced by this extrapolation, as evaluated from the MC simulation, is 
negligible. 
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7.2 Systematic uncertainties 

The systematic uncertainties on the measured cross sections were determined by varying 
the selection cuts and by modifying the analysis procedure. 

For the low-Q 2 sample, the main contribution arises from the uncertainty of ±1 mm in 
the position of the BPC, leading to a ±10% uncertainty in the cross section. 

For the high-Q 2 sample, the dominant sources of uncertainty are as follows. The corre- 
sponding average uncertainties are given in parentheses: 

• the trigger efficiency (IjJ^o); 

• the fiducial volume cut on the electron position was changed by ±1 cm; the SRTD 
alignment was changed by ±2 mm along the Y axis (^35%). A maximum change of 
— 11% was observed in the lowest Q 2 bin; 

• the mass window used for signal extraction was extended by 0.1 GeV (±1.7%); 

• the normalisation of the QED background was changed by ±10% (±2.4%); the maxi- 
mum effect of ±5% was found for the lowest t bin. 

The uncertainty due to the subtraction of proton-dissociative background has been dis- 
cussed in Section 16.21 Additional contributions come from the uncertainties on the in- 
tegrated luminosity, ±2.25%, and on the branching ratios, ±1.7%. Uncertainties from 
the minimum energy requirement of the scattered electron (±0.7%), the elasticity cut 
the selection of the electron and muon samples (±1.2%) and the dependence on 
the MC parametrisations (±0.7%) were also estimated. The total systematic uncertainty 
was determined by adding the individual contributions in quadrature. The correlated and 
uncorrelated systematic uncertainties were evaluated separately and were +g% and ^4%, 
respectively. 

7.3 Dependence on W and Q 2 

The cross section cr^ ^ J ^ p , measured as a function of W and Q 2 for \t\ < 1 GeV 2 , is 
given in Tables Q and El The same cross section, extrapolated to the full t range, is shown 
in Fig. El together with the HI [5] measurements 3 as well as the ZEUS measurement of 
exclusive J/ip photoproduction [9]. The HI measurements are systematically lower than 
the ZEUS data. 

3 In Fig. Et, the HI cross sections, measured at Q 2 values of 3.5, 10.1 and 33.6 GeV 2 , have been rescaled 
to the Q 2 values of 3.1, 6.8 and 16 GeV 2 using the Q 2 dependence of the data measured by HI. The 
systematic uncertainties due to this extrapolation were negligible. 
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The functional form a oc W s was fitted to the ZEUS data; the results of the fit are 
shown in Fig. and in Table 0] No significant variation of 5 with Q 2 is seen. The mean 
value of 5 is 0.73 ± 0.11(stat.)lQ;Qg(syst.). It is consistent with the values found for J/ip 
photoproduction [1] and for p electroproduction at high Q 2 [2]. 

The function a = a ■ (M) li J(Q 2 + Mj^)) n , fitted to the ZEUS data including the 
photoproduction point, is shown in Fig. &p. The resulting parameters are o"o = 77 ± 3 nb 
and n = 2.44 ± 0.08, with % 2 /ndf = 4.1/7. The fit, which takes both the statistical and 
uncorrelated systematic uncertainties into account, describes the data well over the full 
Q 2 range. 

7.3.1 Comparison to model predictions 

Models based on QCD are able to describe exclusive vector meson production at HERA. 
In such models, in the frame where the proton is at rest, the photon emitted from the 
electron fluctuates into a qq state, this qq pair subsequently interacts with the proton 
through the exchange of gluons in a colour-singlet configuration and eventually forms 
a bound meson state. The transverse size of the qq pair depends on Q 2 and on the 
quark mass; for Q 2 > O(10) GeV 2 or for heavy quarks, it is assumed to be considerably 
smaller than the size of the proton. At such distances, the QCD coupling is small and 
perturbation theory can be applied. The QCD factorisation theorem for hard exclusive 
electroproduction of mesons [10] predicts that, in the limit of large Q 2 and fixed x, the 
cross section can be estimated from a hard interaction part calculable in pQCD, the qq 
wave function of the meson and the generalised parton distributions (GPDs) [11] which 
contain information about the correlations of the partons inside the proton and their 
momentum distribution. A rapid rise in the cross section with W is predicted which is 
related to the fast increase of the gluon density inside the proton at small values of x. A 
selection of the available models is compared to the data and discussed below. A more 
complete discussion on the avalaible models is given elsewhere [12]. 

Frankfurt, Koepf and Strikman (FKS) [35] have proposed a model based on the leading- 
order approximation a s \n(Q 2 ). The usual parton distribution functions (PDFs) are used. 
The J ftp wave function is estimated in the non-relativistic approximation. 

In the model of Martin, Ryskin and Teubner (MRT) [36], the calculations are also per- 
formed at the leading order, a s \n(Q 2 ). Assuming parton-hadron duality, the component 
of the cc pair which has the correct spin-parity for the J/ip is used instead of the J/ip 
wave function. The cross section is integrated over the J/ip mass range. The GPDs are 
estimated using the conventional next-to-leading (NLO) gluon distributions. 

Gotsman et al. (GLLMN) [37] have presented a dipole model where the cross section 
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is expressed as the convolution of the wave function of the virtual photon, the dipole 
scattering amplitude and the J/ip wave function. The dipole scattering amplitude is 
estimated at leading order, a s ln(l/x), as the solution of the Balitsky-Kovchegov [38,39] 
evolution equation, including both the linear BFKL terms due to parton splitting and 
nonlinear terms due to recombination of partons in the high-density region at low x. The 
J I xjj wave function is estimated in the non-relativistic approximation. 

The W and Q 2 dependence of the cross sections measured by ZEUS are compared to the 
QCD predictions in Fig.0] As the full NLO corrections have not yet been estimated, all the 
models have significant normalisation uncertainties. Therefore the normalisation was fixed 
using the ZEUS photoproduction data at W = 90 GeV; the different normalisation factors 
are indicated in the figure. The gluon PDFs ZEUS-S [40] for MRT and CTEQ4L [41] for 
FKS were used. The Q 2 dependence of 5 is compared in the insert in Fig. 0Jl All models 
predict a rise of the cross section with increasing W and have a Q 2 dependence similar to 
that of the data. 

7.3.2 Comparison to model predictions for different gluon parametri- 
sations 

The MRT model was used to test three different gluon distributions: MRST02 [42], 
CTEQ6M [43] and ZEUS-S [40], obtained from NLO DGLAP analyses of structure func- 
tion data. In deriving the GPDs from the PDFs, sensitivity to the gluon distribution at 
very low x is introduced. Again, the predictions were normalised to the ZEUS photopro- 
duction measurement at W = 90 GeV. 

Figure El compares the data with the predictions. While CTEQ6M describes the W and 
Q 2 dependence of the data, MRST02 has the wrong shape in W, particularly at low Q 2 . 
ZEUS-S describes the W dependence but falls too quickly with increasing Q 2 . 

The data exhibit a strong sensitivity to the gluon distribution in the proton. However, 
full NLO calculations are needed in order to use these data in global fits to constrain the 
gluon density. 

7.4 Dependence on t 

The differential cross section, da' y * p ^ J ^ p /dt, measured as a function of t in the range 
|i| < 1 GeV 2 , is shown in Table 01 and Fig. EK-d for the high Q 2 sample as well as for 
three Q 2 intervals. A function of the form da/dt = da/dt\ t=0 ■ e~ b ''' was fitted to the 
data and the results of the fit are given in Table |U The slope parameter b is shown 
in Fig. Et as a function of Q 2 and is compared to the ZEUS photoproduction and HI 
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electroproduction values. No significant Q 2 dependence in b is seen over the measured 
range of Q 2 . This behaviour is different from that of exclusive p electroproduction, where 
the b slope strongly decreases with increasing Q 2 , reaching the value of that of the J/ip 
at Q 2 ~ 30GeV 2 [2]. 

In QCD-based models, at high Q 2 , the size of the qq pair in the direction transverse to 
the reaction axis decreases as 1/Q and the t dependence should reach a universal limit, 
independent of the flavour of the quark constituents of the meson [44]. Hence, in this 
limit, the t dependence is given solely by the GPDs of the nucleon. Following this idea, 
the differential cross section was also fitted using an elastic form factor for two-gluon 
exchange, da/dt oc (1 — £/ m 2#)~ 4 5 where m 2 , is the square of the two-gluon invariant 
mass, as suggested by Frankfurt and Strikman [45]. The fit, including both statistical 
and systematic uncertainties, yields m 2 g = 0.55 ± 0.02 GeV 2 and is shown in Fig. 

7.5 Pomeron trajectory 

Soft diffractive processes are described by Regge phenomenology [46] in terms of the 
exchange of a Pomeron trajectory. In hard interactions, where Regge phenomenology 
may not be applicable, an effective Pomeron trajectory may nevertheless be extracted. 
The high-Q 2 sample was analysed to determine the effective Pomeron trajectory. In the 
Regge formalism, the differential cross section can be expressed as 

da/dt oc W iM) ~ l \ (2) 

where the trajectory ajp is usually parametrised as 

OLp{t) = a>]p(0) + a'jpt. 

The effective Pomeron trajectory was determined by fitting Eq. (J2J) to the differential cross 
sections at different t values. The fit was performed in four t bins at Q 2 = 6.8 GeV 2 . Since 
the proton-dissociative process has the same W dependence as the exclusive process, the 
extraction of ajp is not sensitive to this background contribution, which populates the 
high-t region. Therefore the analysis was extended up to |t| = 2 GeV 2 . The fit results are 
shown in Fig. [7| and in Table 03 The parameters of the trajectory, determined from the 
linear fit are: 

a P (0) = 1.20±0.03(stat.)+°;^(syst.); 
a'p = 0.07 ± 0.05(stat.)l^(syst.) GeV~ 2 . 

These values are in good agreement with the ZEUS results from J/ ip photoproduction [9] 
which are also shown in Fig. [7| They are also in agreement with expectations of pQCD- 
based models [47,48], but are not consistent with the trajectory measured in soft diffractive 
processes, ajp = 1.08 + 0.25 t [49,50]. 
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7.6 Decay angular distributions 



The study of the angular distributions of the decay of the J/ ip provides information about 
the photon and J/ip polarisation states. In the helicity frame [51], the production and 
decay of the J/ip can be described in terms of three angles: &h, the angle between the 
J/ip production plane and the lepton scattering plane; Oh, the polar angle, and 4>h, the 
azimuthal angle of the positively charged decay lepton. Under the assumption of SCHC, 
the normalised angular distribution depends only on two angles, Oh and iph = 4>h — 
and can be expressed in the form 

1 dN ^[l+r 4 + (l-3r 4 )cos 2 ^], (3) 



NdcosO 



1 dN 1 r , , 

The spin-density matrix element r°Q represents the probability that the J/ip is produced 
in the helicity-0 state from a virtual photon of helicity or 1. The spin-density matrix 
element r\_ Y gives the probability for the J/ip to be produced in the helicity-1 state from 
a virtual photon of helicity 1 or —1. Assuming SCHC and natural spin-parity exchange 
(NPE) [51], the matrix elements and r\_ l are related by 

^-i = ^(l--o°o 4 )- (5) 



The cross sections at W = 90 GeV are shown in Fig. IHk-f for three intervals of Q 2 . 
Equations (JSJ) and (j3J) were fitted to the data. The values of the spin-density matrix 
elements and rj_ 1: determined from the fits, are given in Table El The measured 
values of r|_ x are consistent with those obtained from Eq. (J3J), also shown in Tabled 
supporting the SCHC and NPE hypotheses. 

Figures 01 and Uni show the cross sections in bins of W and t, respectively. They are quoted 
at the reference value Q 2 = 6.8 GeV 2 . Equation (jHJ) was fitted to the data. The values of 
Tqq, given in Tables [7| and |H1 are consistent with no W or t dependence. 



7.6.1 Longitudinal and transverse cross sections 

The ratio of the longitudinal to transverse cross section, R = cr L /cr T , was calculated as a 
function of Q 2 , W and t from r°o according to the relation 



R 



1 r 04 
1 '00 

e 1 '00 
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which is valid under the assumption of SCHC. 

The values of R as a function of Q 2 are given in Table El and compared with the HI 
results [5] in Fig.[Sfe. The expression R = ((Q 2 /M 2 ^) was fitted to the ZEUS data yielding 
£ = 0.52 ± 0.16(stat.). In QCD-based models, the scale that controls the transverse size 
of the qq fluctuation of the photon may behave differently for and ot- However, in the 
MRT model, <jl and ot have the same W dependence, dictated by the gluon distribution. 
Therefore the ratio is constant. This model correctly describes the rising behaviour of R 
with Q 2 whereas the GLLMN prediction somewhat overestimates it. 

The values of R as a function of W and t are given in Tables and |H1 and shown in Figs. 
andHT%. respectively. 

8 Summary 

The exclusive electroproduction of J/ip mesons, ep — > eJ/tjjp, has been measured with 
the ZEUS detector at HERA for photon virtualities in the ranges 0.15 < Q 2 < 0.8 GeV 2 
and 2 < Q 2 < 100 GeV 2 , for photon-proton centre-of-mass energies in the range 30 < 
W < 220 GeV and for four-momentum-transfer squared in the range \t\ < 1 GeV 2 . 

The cross section of the process 7* p —> J/ip p rises with W as a oc IV s , with a slope 
parameter 5 of about 0.7. This parameter does not change significantly with Q 2 and is 
consistent with that observed in J ftp photoproduction. 

The cross section at W = 90 GeV and over the whole Q 2 range is described by the function 
a oc (Q 2 + M 2 /ip )- n , with n = 2.44 ± 0.08. 

The t distribution, measured for \t\ < 1 GeV 2 , is well described by an exponential de- 
pendence over the range 2 < Q 2 < 100 GeV 2 . The slope parameter, b, is consistent with 
being constant in this range. The mean value is b = 4.72±0.15(stat.) ±0.12(syst.) GeV -2 , 
consistent with that observed in J/ip photoproduction. 

An analysis of the cross sections in the framework of Regge phenomenology yields an 
effective Pomeron trajectory consistent with that measured in J/ip photoproduction. 

The spin-density matrix elements r\_ 1 and are consistent with s-channel-helicity con- 
servation. The ratio of the cross sections for longitudinally and transversely polarised 
photons, R, increases with Q 2 , but is independent of W and t, within the measured 
range. 

The J /ip electroproduction data can be qualitatively described within the framework 
of pQCD that successfully describes J/ip photoproduction data. The data exhibit a 
strong sensitivity to the gluon distribution in the proton. Full next-to-leading-order QCD 
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calculations would allow these data to be used in global QCD fits to constrain the gluon 
density function in the proton. 
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Table 1: The cross sections for the reaction ^y*p — > J/ipp measured as a function 
of W in bins of Q 2 and for \t\ < 1 GeV 2 : (W) and (Q 2 ) are the mean values in 
the indicated ranges; N ee and are the number of events in the signal region 
after non-resonant background subtraction of the electron and muon pairs, respec- 
tively; A ee and are the corresponding acceptances. The first uncertainty of the 
cross sections is statistical and the second systematic. An overall normalisation 
uncertainty of ^1° was not included. 
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Table 2: The cross sections for the reaction ^y*p — > J/ip p measured as a function of Q 2 , for a mean value 
W = 90 GeV and for \t\ < 1 GeV 2 : (Q 2 ) indicates the mean value in the Q range considered; N ee 
and iV MM are the numbers of events in the signal region after non-resonant background subtraction of 
the electron and muon pairs, respectively; A ee and A^ are the corresponding acceptances. The first 
uncertainty of the cross sections is statistical and the second systematic. An overall normalisation 
uncertainty of was not included. 
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Table 3: The differential cross sections for the reaction ^y*p — > J/iftp measured as 
a function oft in bins of Q 2 for a mean value W = 90 GeV . The first uncertainty 
is statistical and the second systematic. 
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(a oc W 5 ) 


b (f oc e" fo l*l) 
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dt lt=0 
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0.15 - 0.8 


0.4 


0.87±0.22lg;{£ 
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5 - 10 
10 -100 


3.1 
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0.65 ±0.171*8 
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1.12 ±0.20lJu| 


4.85 ± 0.24l°:?g 
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2 -100 


6.8 


0.73 ±0.111™ 


4.72 ±0.15lg:ia 
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Table 4: The parameters 5, b and %|t=o measured as a function of Q 2 in the 
range 30 < W < 220 GeV and 45 < W < 160 GeV for the electron and muon 
channels, respectively, and \t\ < lGeV 2 . The first uncertainty is statistical and the 
second systematic. 
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1*1 

(GeV 2 ) 
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0.0 - 0.1 
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0.9 - 2.0 


0.046 
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0.483 
1.123 
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Table 5: The Pomeron trajectory ajp(t) measured in four t bins, in the range 
2 < Q 2 < 100 GeV 2 at a mean value (Q 2 ) = 6.8 GeV 2 . The first uncertainty is 
statistical and the second systematic. 
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Table 6: The spin-density matrix elements, Tqq and r\_ 1} the ratio of cross 
sections of longitudinally and transversely polarised photons, R, and the quantity 
r\_ x — | (1 — r^) measured in bins of Q 2 . The first uncertainty is statistical and 
the second systematic. 
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Table 7: The spin density matrix element and the ratio of cross sections of 
longitudinally and transversely polarised photons, R, measured in bins ofW, in the 
range 2 < Q 2 < 100 GeV 2 at a mean value (Q 2 ) = 6.8 GeV 2 . The first uncertainty 
is statistical and the second systematic. 
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r 04 
'00 


R = o L jo T 


0.0 - 0.1 
0.1 - 0.2 
0.2 - 0.4 
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0.24±0.1ll^ 2 
0.36±0.13l°;°? 
0.19±0.10ig;S^ 
0.16±0.10^o5 


0.31±0.19±g;fg 
0.56±0.30lg;g 
0.23=h0.15±gi| 
0.19±0.14lg;g| 



Table 8: The spin density matrix element and the ratio of cross sections of 
longitudinally and transversely polarised photons, R, measured in bins of \t\, in the 
range 2 < Q 2 < 100 GeV 2 at a mean value (Q 2 ) = 6.8 GeV 2 . The first uncertainty 
is statistical and the second systematic. 
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ZEUS 




Figure 1: The distribution of the events in the muon and electron channels in 
the kinematic plane of Bjorken-x and Q 2 . The events reconstructed using one and 
two measured tracks are shown separately. 
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Figure 2: Invariant mass distributions of the lepton pairs for (a) the low-Q 2 
sample and (b)-(e) the high-Q 2 sample. The shaded histograms are the QED MC 
distributions and the open histograms the sum of the J / ip and QED MC events. 
The small excess of data at low mass is due to background from pions. The error 
bars indicate the statistical uncertainties. 
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ZEUS 




Figure 3: Exclusive J/ip electroproduction cross section (a) as a function of 
W for four values of Q 2 and (b) as a function of Q 2 at (W) = 90 GeV. ZEUS 
photoproduction and HI electroproduction cross sections are also shown. The full 
lines are fits to the ZEUS data of the form (a) a oc W 5 ^ and (b) a oc (Q 2 + 
M 2 ^)~ n . The inner error bars represent the statistical uncertainties, the outer 
bars are the statistical and systematic uncertainties added in quadrature. An overall 
normalisation uncertainty of ^gy° was not included. 
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ZEUS 




Figure 4: Exclusive J/ip electroproduction cross section (a) as a function of 
W for four values of Q 2 and (b) as a function of Q 2 at (W) = 90 GeV . ZEUS 
photoproduction results are also shown. The curves represent the predictions of the 
QCD models MKT, FKS and GLLMN (see text) normalised to the ZEUS photopro- 
duction point at (W) = 90 GeV . The insert shows the parameter 5 as a function 
of Q 2 . The inner error bars represent the statistical uncertainties, the outer bars 
are the statistical and systematic uncertainties added in quadrature. An overall 
normalisation uncertainty of ^1° was not included. 
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ZEUS 




Figure 5: Exclusive J/ip electroproduction cross section (a) as a function ofW 
for four values of Q 2 and (b) as a function of Q 2 at (W) = 90 GeV . ZEUS photo- 
production results are also shown. The data are compared to the MKT predictions 
(see text) obtained with different parametrisations of the gluon density and nor- 
malised to the ZEUS photoproduction point at (W) = 90 GeV . The insert shows 
the parameter 5 as a function of Q 2 . The inner error bars represent the statistical 
uncertainties, the outer bars are the statistical and systematic uncertainties added 
in quadrature. An overall normalisation uncertainty of ~^ly was not included. 
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Figure 6: Differential cross-sections da/dt (a) over the entire Q 2 range and 
(b)-(d) for three bins of Q 2 , for 30 < W < 220 GeV and \t\ < 1 GeV 2 . The full 
lines are the results of a fit to the form da/dt = da/dt\ t=0 ■ e~ b '*' and the dashed 
line is the result of a fit using an elastic form factor assuming two-gluon exchange: 
da/dt oc (1— t/m 2 g )~ A . (e) The slope b, as a function of Q 2 , compared to the ZEUS 
photoproduction and HI results. The mean value ofb is indicated by the horizontal 
line. The inner error bars represent the statistical uncertainty, the outer bars the 
statistical and systematic uncertainties added in quadrature. 
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Figure 7: (a)-(d) Differential cross sections da/alt as a function ofW for fixed 
ranges oft; the full lines are fits to W A ( ap( ^~ l \ (e) Pomeron trajectory: the lines 
are linear fits to the data at (Q 2 ) = 6.8 GeV 2 (full) and to the J/i/j photoproduction 
results (dashed); the dotted line is the soft Pomeron trajectory [52]. The inner 
error bars represent the statistical uncertainty, the outer bars the statistical and 
systematic uncertainties added in quadrature. 
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Figure 8: (a)-(f) Distributions of cos 9 h and ^ in three Q 2 bins; the curves are 
the fits to Eqs. ([3J) and Q). (g) Ratio R = olI®t as a junction ofQ 2 ; the full curve 
is the result of the fit to the ZEUS data while the dashed and dotted curves are the 
predictions of the MRT and GLLMN models, respectively. The inner error bars 
represent the statistical uncertainty, the outer bars the statistical and systematic 
uncertainties added in quadrature. 
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Figure 9: (a)-(e) Distributions of cos 8h for five bins ofW; the curves are the fits 
to Eq. ([3J). (f) Ratio R = o^jo^ as a function of W; the dashed line is the MRT 
prediction and the full line is the result of a one-parameter fit. The error bars are 
statistical (inner) and total (outer). 
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Figure 10: (a)-(e) Distributions of cos 9 h for four bins oft; the curves are the 
fits to Eq. (f) Ratio R = o~lIo~t (is a function oft; the full line is the result of 
a one-parameter fit. The error bars are statistical (inner) and total (outer). 
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